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^ ■ Abstract 



We study the phase diagram of the Ising-Kondo lattice with transverse mag- 



o 

0> , netic field as a possible model for the weak-moment heavy-fermion compound 



URu2Si2, in terms of two low- lying / singlets in which the uranium mo- 
ment is coupled by on-site exchange to the conduction electron spins. In the 



c . 

O . mean-field approximation for an extended range of parameters, we show that 

O 

the conduction electron magnetization responds logarithmically to /-moment 
formation, that the ordered moment in the antiferromagnetic state is anoma- 
lously small, and that the Neel temperature is of the order observed. The 
model gives a qualitatively correct temperature-dependence, but not magni- 
tude, of the specific heat. The majority of the specific heat jump at the Neel 
temperature arises from the formation of a spin gap in the conduction elec- 
tron spectrum. We also discuss the single-impurity version of the model and 
speculate on ways to increase the specific heat coefficient. In the limits of 
small bandwidth and of small Ising-Kondo coupling, we find that the model 
corresponds to anisotropic Heisenberg and Hubbard models respectively. 
PACS numbers: 71.10.Fd, 71.27.-ha, 71.70.Ch, 75.30.Mb 



I. INTRODUCTION 



As the large specific heat jump of URu2Si2 at its Neel temperature of 
Tat = 17.5 K is difficult to reconcile with its anomalously small staggered 
magnetizationBI, some authorJl"! have speculated that there is some hid- 
den order parameter. Recent neutron experiments indicate, however, that 
the weak Bragg peaks of the ordered phase break time reversal symmetry as 
would be the case if the order were magnetic dipolaiH. In addition the hidden 
multispin order parameters proposed to account for the large specific heat 
jump have been shown to give the opposite field-dependence for the Bragg 
peaks to that observed in experimentfl. To date, neither the small moments 
(about 1% of the full uranium momenlSi) nor the large specific heat jump! 
AC/C ~ 1 at T/v have been explained. In this paper, we present a model 
which may account for the small moments, and discuss additional require- 
ments needed to attain sufficiently large specific heat jump AC and linear 
specific heat coefficient 7. 

The sharp, dispersive spin excitations in URu2Si2 are longitudinal, and 
are well described in terms of a transition between two crystal-field singletsSJl; 
no sharp higher energy levels have been observed, consistent with the highly 
anisotropic susceptibility. The low transverse susceptibility (x^/x^'^ ~ 5 as 
T ^ 0) indicates that the next higher degenerate crystal-field levels are at 
energies an order of magnitude higher than the A = 2.4 THz of the first 
excited singlet state. A broad longitudinal continuum has been observed at 
higher energies, but we focus here on the sharp crystal-field-like features of the 
spectrum. Although a doublet crystal-field ground state has been proposed!! 
for U diluted in ThRu2Si2 with its different cell parameters, we follow the 
evidence from neutron scattering that URu2Si2 has two low- lying singlets 
separated by a gap A. 
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In an attempt to describe the magnetism of URu2Si2, we consider the 
following Hamiltonian which is an Ising-Kondo lattice model with transverse 
field: 

k « « 

where ip\ = (cj-|-,ctj^) (with as its Fourier transform), and creates a 
conduction electron of spin a =t,| at lattice site i] is the conduction 
band energy, with width 2W . The 5=1/2 operators Si, normalized so 
Sf = 3/4, correspond to the two-level system consisting of the two low- lying 
singlets which are separated by an energy A, and we thus regard it as the 
spin operator for / electrons: the isotropic Kondo exchange / Si ■ {ip} ^^il^i ) 
reduces to the above Ising form when the true /-spin is projected onto the 
lowest two singlet states. Details about the crystal-field splitting are given in 
Appendix ^. 

In Ref. [|lO|, parameters of this Hamiltonian were estimated that agreed 
with the expected 6d conduction-electron bandwidth 2W ~ 500 THz, on-site 
exchange / ~ 25 THz comparable to that in other heavy-fermion systems, 
and the singlet-singlet gap A ~ 2.4 THz observed in neutron scattering at 
low T. These authors also reported an order-of-magnitude mass enhancement 



from spin fluctuations above T^, but as discussed in Sec. |V|, a lowest order 
perturbative calculation shows that the mass enhancement is small unless 
I/W is 0(1). 

While URu2Si2 does not necessarily correspond exactly to half-filling, in 
this paper we restrict attention to the half-filled case, to one-dimensional 
(ID) and centred tetragonal0 (ct) lattices, and a conduction-electron band 
constructed from inter-sublattice hopping: 



E ^fc^K = E 4s- + h-c. 



where {ij) denotes bonds between nearest neighbours on opposite sublattices. 
Specifically, 

-2tcos{ka), ID 
—8tcos{kxa/2) cos{kya/2) cos{kzc/2), CT 
so that the bandwidth 2W is At (16t) in the ID (ct) lattice. Here the conven- 
tional CT unit cell, considering the CT lattice to be BCT, measures a x a x c. 

This paper proceeds as follows. We study the model mainly in the mean- 
field approximation, and find an extensive low-moment regime as well as rea- 
sonable Neel temperature, although the specific heat, while showing a jump 
at Tn, is not of sufficient magnitude. Using a single- impurity model for the 
specific heat, wc speculate on ways to obtain a large specific heat coefficient 
7. We also find that in the limits of small bandwidth and low coupling, direct 
comparisons with other electron correlation models are possible. 



II. MEAN-FIELD APPROXIMATION 



We present mean-field results for our Hamiltonian which naturally provide 
a region of parameter space in which the required small moments occur. 

Taking as our mean-field ansatz Sj = {—lynif + SSj for the /-electron 
spins and 'ipj^a^ipj = s| = {~iy~^^mc + Ssj for the conduction-electron 
spins, where (—1)-^ = exp(zQ • Rj) alternates between opposite sublattices, 
our Hamiltonian becomes, upon dropping the product of fluctuations. 



/ 
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where M is the number of uranium atoms, and k runs over the magnetic 
Brillouin zone. Minimizing the free energy, we obtain self-consistent mean- 
field equations for the /- and conduction-electron moments irif and rric, given 



here for bands satisfying the generahzed nesting condition = — e^_,_Q 
for all k (which occurs, for example, for nearest-neighbour inter-sublattice 
hopping in ID, square, and CT lattices): 

nif = ^^tanh(/?A/2) , 
_ Inif y-^ tanh(/3e^/2) 



where ite^ = ibwe?.-|- Q^IrajY and A = + {Irric)'^ are the mean-field 



energy levels, and fi is the inverse temperature. We note that the mean- 
field moments ruf and rric both have mean- field critical exponent ^ and are 
therefore proportional to each other near the Neel temperature T/v, which is 
obtained by the solution of 

.u/A/oi.^^ /2 tanh(e^/2fcB7V) 
The specific heat per uranium atom is 



k ^ 



N 

k 
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„ /2 dm?f 
^' kl8T dT 



2A;|T dT 



sech2(/3e^V2) 
sech2(/3A/2) . 



For comparison to experiment, we note that a staggered magnetization 
nif = 0.5 in our Sf = ^ model corresponds to a true moment of 1.2/xb, as 
this is the observed matrix element of the true operator between the two 
smg letH. In Fig. H we show that small /-moments naturally arise from the 
mean-field approximation in a region of parameter space not too far from the 
values of Ref. quoted above. Mean-field theory also gives a Neel temper- 
ature of appropriate order. We expect that in a full calculation fluctuation 
effects would reduce both the moments and Neel temperature. Thus the val- 
ues of the parameters could differ from those derived by matching mean-field 
theory to experiment. 



It is difficult within mean-field theory to obtain both the Neel temperature 
and moments of experiment since we do not have freedom to simply scale all 
parameters: the singlet-singlet spacing cannot be very different from A = 
2.4 THz. 

Since the Fermi surface is perfectly nested (in the generalized sense de- 
fined above) by the antiferromagnetic wave-vector — this is true for nearest- 
neighbour inter-sublattice hopping in ID, square, and CT lattices — we find 
that, as expected, the conduction electrons exhibit logarithmic response to the 
field produced by the /-system: if Imj <^ W, then rric ~ m/| logm/|'^ in d 
dimensions. This behaviour is caused by the sharp corners found in the Fermi 
surface only at half-filling, and is therefore not present away from half-filling, 
or for more general Fermi surfaces. 

The temperature-dependence of the /-moment and specific heat for se- 
lected parameters are given in Fig. ^ for the ID case, showing a jump 
AC/C ~ 15% with C too small by a factor of 6 to fit experiment. The 
jump arises mainly from the formation of a spin gap Imj in the conduction 
electron spectrum. 



III. ANALYTIC BEHAVIOUR IN TWO LIMITS 

We now study our Hamiltonian in two relatively simple limits, namely 
those of small bandwidth and small coupling which to first order lead, respec- 
tively, to the Heisenberg model with Ising anistropy and the Hubbard model. 
It is clear that the inclusion of additional terms in our Hamiltonian to model 
the real system will dominate some of the terms we find in the effective Hamil- 
tonians in these limits. For example, had we included a nearest-neighbour 
Coulomb interaction between conduction electrons, the U term we obtain for 
the Hubbard model in Sec. Ill B below would be neg hgiblell. So while we do 
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not believe that these hmits pertain to the real material URu2Si2, we study 
them to obtain a more global picture of our Hamiltonian and hopefully under- 
stand more about its behaviour in the real, relatively complicated, parameter 
regime. Fig. ^ depicts a summary of the results of this section. 

A. Small Bandwidth Limit 

If t = 0, then our Hilbert space breaks up into an independent-site descrip- 
tion which is easily diagonalized to give the degenerate energy-level diagram 
for each site shown in Fig. ^ In the small t limit, we are restricted to the 
ground state doublet in Fig. |^, in the sector of the Hilbert space in which 
there is exactly one conduction electron per site. Taking the standard strong- 
coupling approach for 



t«iv//2 + 4A2-lA, 

we obtain in second-order degenerate perturbation theory (similar to that 
for the large-[/ Hubbard modelEl) the Heisenberg Hamiltonian with Ising 
anisotropy: 

where 

'' = p(pTTaw(''^""'^'^''^') 

In the limit A// ^ we have J^'^ — > 0, which gives the pure Ising Hamil- 
tonian. The opposite limit //A — > results in the Heisenberg Hamiltonian 
with 

32t2A 



and Ising anisotropy 



A3 



B. Low Coupling Limit 

Using a path-integral method, we integrate out the localized /-spins in 
the limit / ^ A, and we find the imaginary-time effective Lagrangian for the 
conduction electrons to be 

= E 4. (^^ + ^.^) ^a:. - ^ E + E {d.s^^ir)f + . . . 

where • • • represents terms of higher order in derivatives and in //A. 

With just the first correction term (the one proportional to /^/A), noting 
that = \nj — ^nj|njj, we have an effective Hubbard Hamiltonian with 

8A 

which in the strong-coupling limit t <^ U gives a t-J model with 

4i2 _ 32t2A 

which agrees with our small t limit in the special case of I/A shown in 
Eq. (H). 

At half-filling, the Hubbard model is believed, on the basis of mean-field 
theory, to exhibit a Mott-Hubbard transition at infinitesimal U to an antifer- 
romagnetic phase0. In d dimensions, the staggered magnetization is 



rric oc — exp 



-{at/U) 



l/d 



with the singularity arising from the Fermi surface corneral^. For the model 
of current interest in the {d = 3) CT lattice, we determine the coefficient in 
the exponential to be a = 3tt'^/2. 
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We now interpret 

drs'j = [H, s'j] ^ [Ho, s'j] + 0(//A), 

where Hq is the Hamiltonian with 7 = 0. This results in an effective Hamil- 
tonian which has interaction terms 

j2 ^2j2 ( / 

j I (ij) ^ " 

where J2{ijk) = J2j(i=j-i,k=j+i) (vahd only in ID) and the hermitian conju- 
gate h.c. applies to all terms to the left within the brackets. 

We note that in the small t limit, in which we project onto the singly- 
occupied subspace, all but the first two terms disappear and the result agrees 
precisely with the Ising anisotropy found in Eq. (U). 

IV. LARGE EFFECTIVE MASS 

While the mean-field approximation gives experimentally interesting (and 
hopefully accurate) values for the zero-temperature moment and Neel tem- 
perature, in addition to not predicting a large specific heat jump, it does 
not account for another fundamental property of URu2Si2: the large value 
of 7 180 mJ/mol-K^, the zero-temperature intercept of C/T above Tjy, 
i.e. the heavy- fermion mass. This is because, once the moments have gone 
to zero, the mean-field specific heat simply becomes that of free conduction 
electrons plus a Schottky term from the /-moments. 

One possible way of improving on mean-field theory above T^r is to ex- 
trapolate known or conjectured results on the single-moment version of the 
problem to the lattice version. That is, we consider the Hamiltonian of Eq. (|l]) 
with a single spin, located at the origin. This model may be recognized as a 
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particular realization of the well-studied problem of a spin coupled to a heat 
bath. Assuming a spherically symmetric dispersion relation, and expanding 
the fermion fields in spherical harmonics, only one harmonic interacts with 
the impurity spin so the problem is effectively one-dimensionalEl. This one- 
dimensional spin-fermion problem may be bosonized at low energies, giving 
the spin-boson problem. This model has been used to study the effect of 
dissipation on tunneling. The dimensionless strength of the dissipation is 
measured by I/W; A corresponds to the tunneling matrix element. Since we 
are apparently in the weak dissipation limit, I/W <^ 1, the spin has a unique 
ground state with (5"^) ^ 0. (At stronger dissipation there are two ground 
states with (S^) = itm ^ 0.) The spin-boson problem in turn is equivalent, at 
low energies, to the Kondo problem. The weak dissipation case corresponds 
to antiferromagnetic Kondo coupling with a screened ground state. 

According to Ref. this connection with the spin dissipation problem 
gives (57 oc f{I/W)/A with f{x) as x ^ and /(I) 1. Using 

perturbation theory to second order in I/W, we indeed find that 

where po = 1/2W is the density of states at the Fermi level. But we need 
7 ~ to fit experiment; with I/W ~ .1 this is two orders of magnitude 
too small. In order to obtain a large enough 7 we would apparently have to 
choose I/W ~ 1, in which case the moment would not be small according to 
Fig. |1|. 

Thus, making the non-interacting impurity approximation, we cannot ex- 
plain the large 7 value without assuming a large value of I/W. It is possible 
that treating our model more accurately, including rkky interactions, will 
lead to a sufficiently large 7. Alternatively, the model may be missing some 
important physics. Recall that we have thrown away the uranium crystal field 

10 



levels which are an order of magnitude higher; reinstating them brings the 
spin-flip part of the Kondo interaction back into play. 

The huge mass enhancement of the charge carriers in many heavy- fermion 
compounds is often explained in terms of the Kondo effect. Besides the bare 
conduction-electron bandwidth, the Kondo effect introduces a much smaller 
energy scale, namely the Kondo temperature Tk, below which the local mo- 
ment degrees of freedom are frozen. The same small energy scale also gives 
rise to an effective narrow band with an enhanced density of states: if the 
Kondo temperature is of order A (the crystal-field splitting) or larger, the 
spin-flip parts of the Kondo coupling can renormalize to the strong-coupling 
fixed point producing a 7 per ion of 0{1/Tk)- In URu2Si2, we expect this 
renormalization to be cut off at a scale of order A; however, one still ex- 
pects the same kind of Kondo screening process over a wide range of energy 
scales from the order of the bare bandwidth down to 0(A). While this effect 
is clearly not included when we write down the low-energy effective Ising- 
Kondo lattice model, a semi-phenomenological way of treating it is to use our 
model with a greatly reduced effective bandwidth of 0{Tk)- This effective 
low-energy theory is valid at low energy scales after the higher crystal-field 
levels have been integrated out and the associated Kondo effect has been taken 
into account. 



V. DISCUSSION 

We have shown that within the mean-field approximation, small moments 
are predicted for a range of parameters because of a logarithmic response 
of the conduction electrons at half-filling. A specific heat jump is obtained 
mainly from the formation of a conduction-electron spin gap at Tjv. 

The question remains as to the nature of the dynamical narrowing of the 
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bandwidth. 

We have given comparisons of our model in various hmits to other models, 
such as the Hubbard, Heisenberg, and Ising models. Because anisotropy is 
maintained in these limits, we expect that the one-dimensional system will 
still order at zero temperature. 
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APPENDIX A: DETAILS OF THE CRYSTAL-FIELD SPLITTING IN URugSis 

We take the uranium ion to have total angular momentum J = 4, and the 
crystal-field Hamiltonian which splits the 9-fold degeneracy has the form0 

HcF = Biol + Biol + Blot + Biol + Btot (Al) 

where are the Steven's angular momentum operators for centred tetrag- 
onal symmetry; B"!^ are the coefficients for URu2Si2. With the coefficients 
we choose in an attempt to match the neutron-scattering experiments, this 
Hamiltonian results in two low-lying singlets as well as two high-lying doublets 
and three singlets. We represent what we consider the important structure in 
Fig. |5|. The two low-lying singlets are 

|0) = |r,\) ^ 6(|4, 4) + |4, -4)) + 7|4, 0), (A2) 
|l) = |r,2)^^(|4,4)-|4,-4)) (A3) 
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where |J, mj) are the states of total angular momentum J and azimuthal 
quantum number mj. The operators bring |0) and |1) into the high-energy 
states, while (0|J^|1) = 8e/\/2- Writing s for the conduction electron spin, 
the Kondo interaction J ■ s = J^s^ + ^{J^s^ + J^s+) reduces to just its Ising 
part J^s^ when considering only the two low- lying singlets. Now we define a 
dimensionless and normalized spin-^ operator such that = (8e/^/2)5'^. 
In order to implement the separation A between the two low-lying singlets, 
we introduce an operator which is diagonal in this subspace: 

S^\0) = -l\0); 

Thus we obtain the terms ^f-^i + Sf in our Hamiltonian (|l|). 
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FIGURES 
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FIG. 1. Zero-temperature /-moments and Neel temperatures for ID (dashed lines) and CT 
(solid lines) lattices as functions of the bandwidth, singlet spacing, and Ising-Kondo coupling. In 
each case we fix two of the parameters and vary the other. Experimental values m/ = 0.04/xb and 
T/v = 17.5 K are indicated on the axes. 
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FIG. 3. Since overall scaling is irrelevant, the space of points {t, I, A) is two-dimensional. Here 
the points are represented in such a way that t : I : A = diA '■ dtA '■ dti where dab is the distance 
from the point to the line joining a and b. The regions we have examined are the shaded ones. 
(For the purposes of this figure, we interpret a <C 6 as a < b/3.) 
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FIG. 4. Spectrum of the independent-site Hamiltonian for t = 0. Here the first 

kct represents the conduction electrons: \a) = cj.|0) for a =],[', \ U) = cjcj |0); 
|0) is the "vacuum". The second ket represents the state of the locahzed / electron: 

I y) = \ve); \ \) = \vn-e); \ = \Vn/2); \ ^) = \v-n/2); \ \) = \v-e); \ /) = Ive-n), and 

\V24>) = ifjcoscj) + /|sin(/))|0) is the state in which the /-electron spin lying in the x-z plane 
makes an angle 2(f) with the positive z-axis measured toward the positive x-axis — note that 
(r/20|772</)') = cos(^ — <f)'). It is easy to see that here 6 = tan~^(2A//) is the tilt produced by the 
"transverse magnetic field" A. 
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FIG. 5. Spectrum of the crystal-field Hamiltonian, showing states related to the low-lying 
singlets by application of the operator J. 
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